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Abstract. During development, morphogen gradients precisely determine the 
position of gene expression boundaries despite the inevitable presence of fluctuations. 
Recent experiments suggest that some morphogen gradients may be interpreted prior to 
reaching steady-state. Theoretical work has predicted that such systems will be more 
robust to embryo-to-embryo fluctuations. By analysing two experimentally motivated 
models of morphogen gradient formation, we investigate the positional precision of 
gene expression boundaries determined by pre-steady-state morphogen gradients in the 
presence of embryo-to-embryo fluctuations, internal biochemical noise and variations 
in the timing of morphogen measurement. Morphogens that are direct transcription 
factors are found to be particularly sensitive to internal noise when interpreted prior 
to steady-state, disadvantaging early measurement, even in the presence of large 
embryo-to-embryo fluctuations. Morphogens interpreted by cell-surface receptors can 
be measured prior to steady-state without significant decrease in positional precision 
provided fluctuations in the timing of measurement are small. Applying our results to 
experiment, we predict that Bicoid, a transcription factor morphogen in Drosophila, is 
unlikely to be interpreted prior to reaching steady-state. We also predict that Activin 
in Xenopus and Nodal in zebrafish, morphogens interpreted by cell-surface receptors, 
can be decoded in pre-steady-state. 
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1. Introduction 

A central paradigm for spatial positioning in embryonic development is the morphogen 
gradient: signaling molecules (morphogens) form a concentration gradient across 
the system with gene expression boundaries determined by whether the morphogen 
concentration is greater or less than a particular threshold [H [2]. Despite the 
apparent simplicity of the morphogen concept there are still many unanswered questions 
including: how morphogens are transported through the embryo [31 HJ El El E] ; how the 
system measures the morphogen concentration [U [9j [10] ; what mechanisms control the 
timing of morphogen interpretation [21 |TTJ [T2J [13] ; and how robust morphogen gradients 
are to fluctuations [TH [151 EE] • 

Experimental and theoretical approaches often assume that morphogen gradients 
are in steady-state prior to measurement [T6l [T71 [18]. However, recent experimental 
evidence suggests that some morphogen gradients may be interpreted before reaching 
steady-state including Sonic Hedgehog (SHH) [101119], Hedgehog |8|, Activin [20], Bone 
Morphogenetic Protein (BMP) [TBI [21], Nodal [22] and, controversially, Bicoid (Bed) 
[23] . Theoretical analysis suggests that such pre-steady-state readout can provide 
reliable positioning of gene boundaries [211 [21] and it may even be preferable to steady- 
state measurement if there are large variations between embryos in the morphogen 
production rate [23]. Furthermore, time- varying morphogen concentrations offer the 
additional advantage that they can define the expression of multiple genes at similar 
spatial positions but at different times [201 IS]. 

Morphogen gradients are known to provide very precise positional information 
about cell fate despite embryos being inherently noisy [25J [261 123 (211 [29] . The gradient 
must therefore be robust to external embryo-to-embryo variations, such as differences 
in size [251 121] and morphogen production rates [TSJ. Internal fluctuations, caused 
by processes such as protein production, diffusion and degradation, are present in all 
biological systems and limit the precision of biochemical signaling [T5l [TBI [30] . To 
reduce internal fluctuations, the embryo time-averages the morphogen gradient [T5], 26J. 
However, in pre-steady-state, the concentration varies systematically as a function of 
time and so variations in the onset of averaging (as well as in the averaging period) could 
reduce precision [23]. For the first time, our work considers all the above fluctuations 
to build a more complete understanding of the precision of pre-steady-state morphogen 
gradient interpretation. 

We analyse two experimentally relevant models of morphogen gradient formation 
[31] : morphogen production restricted to a localised region with global diffusion and 
degradation; and morphogen production by a spatially distributed mRNA gradient [6]. 
In the latter case, morphogen is produced throughout the embryo but is confined to the 
cell/nucleus in which it is produced. We also apply our analysis to a simplified model 
of the Bed gap-gene regulatory network [231 [32] to demonstrate that our approach 
extends to more realistic systems. We find that internal fluctuations present a major 
barrier to precise pre-steady-state interpretation of morphogen gradients in both models. 
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In particular, if the morphogen is a direct transcription factor (with nanometer-scale 
DNA binding sites on the target genes) then internal fluctuations can make pre-steady- 
state decoding especially imprecise. We apply our analysis to Bed, which has recently 
provoked debate over whether it is interpreted prior to steady-state [331 EH- We find 
that a Bed diffusion rate of lpm 2 s' 1 (as suggested by Bergmann et al. [231 E3]) is 
inconsistent with precise position determination within the known time window of Bed 
target gene expression due to error from internal fluctuations. 

Our results can be applied to morphogens interpreted by cell-surface receptors, 
such as Activin in Xenopus and Nodal in zebrafish. The effective measuring measuring 
volume is large for such morphogens (roughly the cell size) and hence the effects of 
internal fluctuations on the precision with which the morphogen is interpreted by the 
cell are small. Provided that fluctuations in the onset time of measurement are small, we 
predict that pre-steady-state measurement precision for such systems is not significantly 
different from that in steady-state, meaning pre-steady-state interpretation is possible. 

2. Models 

We examine two models of morphogen gradient formation: localised morphogen 
synthesis with long-range diffusion and spatially uniform linear degradation (SDD); 
and spatially-distributed mR NA gradients (MRG) with no extra-cellular/extra-nuclear 
diffusion. The SDD model is the prevailing model for morphogen propagation [35] and 
is representative of diffusion-dominated systems; for example, the Activin gradient in 
Xenopus forms through passive diffusion [36J. However, recent experiments on Bed in 
Drosophila, suggest that the bed mRNA itself may form a spatially distributed gradient 
P, EI]. The MRG model provides a simple description of such gradient formation. 
Our models encapsulate the properties of these mechanisms for morphogen gradient 
formation, enabling a more general discussion of the (dis) advantages of pre-steady-state 
interpretation. 

SDD Model: In steady-state the SDD model results in an exponentially decaying 
concentration profile. Bed and Activin are observed to have such concentration profiles 
in steady-state [251 [37]. The morphogen concentration p sdd (~x~,t) at position x and time 
t is described by 

d t p sdd (x, t) = DX7 2 p sdd (x, t) - //p««(x, t) , (1) 

where D is the diffusion constant and p the rate of morphogen degradation. The 

morphogen is produced along the plane at x = with constant flux J, leading to 

the boundary condition Dd x p(~x, t)\ x =o + J = 0. In the absence of fluctuations, the 
general solution of (JTJ is [23] 

p add (x,t) = {l-f(x,t))p e-*' x «", (2) 

where \ sdd = x/DjJi is the decay length, p = JX sdd /D, f{x,t) = erfc(z_)/2 + 
e 2x / Xsdd erfc(z + ) /2, erfc(z) is the complimentary error function, z T = (2Dt/\ sdd =)= 
x) I y/ADt, and x, X sdd <^L (where L is the system length). In steady-state the morphogen 
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concentration is p e~ x ^ Xsdd . The morphogen concentration depends only on the time after 
initiation and the linear distance from the plane at x — 0. Solving ([1]) using a geometry 
similar to that of the Drosophila embryo results in no significant difference in the 
morphogen profile from the approach described above [23], justifying our methodology 
MRG Model: We assume the underlying mRNA gradient is fully formed before 
measurement begins, consistent with experimental observations [61 [31]. Including 
spatially uniform morphogen degradation, the morphogen concentration p mr9 (x, t) is 
described by ^p mr . 9 (x, t) = j — up mrg , where j = j e~ x / Xmrs is the morphogen production 
rate per unit volume and v the rate of morphogen degradation. In the absence of 
fluctuations, the concentration profile is given by 



where p\ = j jv. The gradient is always exponential with decay length X mrg 
(independent of z/), but with increasing amplitude as time increases. 

Figure 1A shows the evolution of ([2]) and ([3]) with time at fixed position, xt, using 
p = v and po = Pi- Pmrg increases more rapidly initially than p S( id as morphogen is 
produced locally at xt- As morphogen begins to diffuse into the measuring region, 
Psdd rapidly increases. Finally, at large times the two models tend to the steady- 
state concentration. As we now demonstrate, the different evolution of the morphogen 
concentrations with time (Figure 1A inset) can have an important impact on the 
positional precision of the two models in the presence of different sources of noise. 

3. Results 

3.1. Internal Biochemical Fluctuations 

Morphogen production, diffusion and degradation are stochastic processes, resulting in 
internal fluctuations in the morphogen concentration about the threshold position. We 
are interested in the positional error, uj int) due to such fluctuations. However, before we 
calculate the spatial precision of the two models in the presence of internal fluctuations, 
we first consider when and how the system interprets the morphogen gradient. At early 
times in one possible scenario, morphogen interpretation could be suppressed by an 
inhibitory protein distributed uniformly across the system. The inhibitor could, for 
example, compete with the morphogen and its cofactor by binding cooperatively to 
morphogen transcription factor binding sites. When the inhibitor density decreases 
(for example, by degradation) below a threshold, repression will be released and 
morphogen interpretation switched on at that particular time throughout the system. 
Even after this switch on, the system cannot make instantaneous measurements of the 
morphogen gradient as the probability of a receptor/DNA binding site being occupied 
by a morphogen molecule is often small [T5 l 126 1 138] . To increase precision, developmental 
systems sample the morphogen gradient over a period of time, referred to here as the 
'time-averaging window' (TAW). The period of the TAW, r m , can be controlled by the 
lifetimes of the downstream products and feedback between the morphogen and these 
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downstream products. For example, chick neural development is partly controlled by a 
gradient of SHH, but cells become gradually desensitized to ongoing SHH exposure due 
to SHH-dependent upregulation of a protein that inhibits SHH signaling |10j . 

During morphogen interpretation, the time between independent measurements 
of the morphogen concentration scales as Ti n d ~ (Ax) 2 / Do, where Ax is the linear 
dimension of the measuring region [151 EO] and D is the local diffusion constant in 
the vicinity of the receptor/DNA binding site. Note that while there is no extra- 
cellular/extra-nuclear diffusion in the MRG model, the morphogen can diffuse within 
the cell/nucleus. 

In our simple model, the statistics of the particle number at a given position due to 
internal fluctuations are Poissonian [T5l 1391 |4"0] . We have confirmed this result in pre- 
steady-state using three-dimensional stochastic simulations for both the SDD and MRG 
models, Figure IB. The stochastic simulations involved simulating the propagation of up 
to 3000 particles using Monte Carlo techniques and measuring the particle occupation 
at specified sites (see [15] for details). Furthermore, we assume that r m <C r, meaning 
that the gradient varies comparatively little over the averaging period, which begins 
at time r. This assumption is reasonable for Bed, where target gene expression begins 
around 90 minutes after initiation but the TAW is restricted by the relatively short time 
between successive nuclear divisions [26]. Hence, we calculate oj int using [T5] 

(Pr,r m (XT)) 1/2 
{d x p T ,r m (x T )) 

where (...) denotes ensemble averaging, p TTm (x) is the average morphogen concentration 
over the period r < t < r m (see Appendix A), and A is a constant determined by 
details of the time-averaging network (see Appendix A). To test the validity of (HI), we 
compared the positional error found from stochastic simulations of the SDD model 
on small systems (see [15] for details) with (j4|). As can be seen from Figure 1C, 
the agreement is good. Four properties of (Jl|) are particularly important: (i) longer 
averaging times result in reduced positional error from internal fluctuations; (ii) low 
concentration results in increased precision due to reduced fluctuations in the density; 
(iii) shallow gradients result in large positional error; (iv) the constant A oc (Ax) -1 / 2 in 
three dimensions, see Appendix A, so smaller measuring volumes result in larger error 
due to internal fluctuations. Therefore, we expect morphogenetic transcription factors 
(Ax nanometre scale, the typical size of a DNA binding site) to have larger error than 
morphogens interpreted by cell-surface receptors (Ax micrometre scale, the typical size 
of a cell). In the latter case, although the receptors themselves are small, the effective 
morphogen measuring volume is comparable to the cell size [301 HI]- 

At early times the morphogen gradient is flat, accentuating any concentration 
fluctuations into large positional error, see Figure 2A (Figures 2 and 3 use parameters 
consistent with the direct transcription factor Bed [231 [261 [35] unless stated otherwise). 
Early measurement in both the SDD and MRG models comes at a major cost to 
precision. The SDD model is particularly imprecise at short times, even including 



u int Pd Ar m 1 — — — T7- , (4) 
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spatial averaging as insufficient morphogen has diffused into the measuring region 
for reliable measurement. The MRG model is more precise than the SDD model at 
early times as morphogen is produced locally at the threshold position, resulting in a 
steeper morphogen profile, see Figure 2A. However, the precision is still worse than 
at steady-state. It is clear from Figure 2A that, given a particular final steady-state 
profile, improved precision is achieved through fast morphogen production, diffusion and 
degradation (i.e. rapid approach to steady-state). Our results are robust for a range 
of biologically reasonable parameters (we explored a range of morphogen production 
rates (J = 0.1 — decay lengths (A = 70 — 120/xm) and threshold positions 

(xt = 150 — 250/im)). In summary, pre-steady-state measurement does not provide 
increased robustness to internal fluctuations and, indeed, it can result in significantly 
increased positional error when the morphogen is a direct transcription factor, due to 
the small measuring volumes involved. 

Rather than the absolute morphogen concentration used above, an embryo could 
alternatively interpret a normalised concentration [U [26]. For example, the nuclear 
Bed concentration appears to be approximately constant between cycles 9-14, although 
the extra-nuclear Bed concentration may vary [35]. In such systems, the numerator 
in (j3J) would be constant (here taken to equal the steady-state value) with changes in 
precision determined by the slope of the morphogen gradient at the threshold position. 
At early times the gradient is relatively flat at typical threshold distances, resulting in 
large positional error. The positional error decreases with time as the gradient becomes 
steeper at the threshold position. We also note that, compared with fluctuations in 
steady-state normalised concentrations, the effects of external and TAW fluctuations 
(see below) are increased by interpreting normalised morphogen concentrations in 
pre-steady-state. Therefore, considering normalised concentrations does not alter our 
conclusions. 

3.2. TAW Fluctuations 

Variations in the TAW affect the spatial precision of pre-steady-state morphogen 
gradients. The effect of different averaging onset times is illustrated for both the SDD 
and MRG models in Figure 2B. We denote the positional error due to such fluctuations 
(St) by ujtaw- Fluctuations in r m can also alter ujtaw (and u) int due to an altered 
averaging period). However, in the regime r m <C r, such fluctuations make little 
difference to our results (data not shown) and are thus neglected. In this regime, the 
positional error due to fluctuations in the averaging onset time at position xt is (for 
St < r) 



where g sdd = Td t f(x T , r) [1 - f(x T , r) + X sdd d x f(x T , r)] 1 and g mrg = vT(e UT - 1) 1 (see 
Appendix B). Interestingly, we see that g mrg is independent of the threshold position 
and of Xmrg, whereas g sdd is sensitive to all the kinetic parameters and the threshold 




(5) 



position. The positional independence of ^taw 



arises because the concentration profile 
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in the MRG model is always exponential, though with a time-dependent amplitude. 
Since morphogens typically control more than one target gene, oj^ r ^ w being independent 
of position could be a possible advantage for the MRG model in pre-steady-state. 

We now investigate the behaviour of (jJtaw numerically (see Appendix A). We first 
analyse how the positional error depends on the magnitude of the TAW fluctuations. 
An important result, clear from Figure 2C, is that the MRG model is typically more 
robust to TAW fluctuations than the SDD model. This result is general if both models 
have similar kinetic parameters, such as the decay length and degradation rate. We see 
that TAW fluctuations can be a significant source of positional error in pre-steady-state, 
even at small values of 5t/t. However, the positional error not only depends on the 
magnitude of the TAW fluctuations, but also on the averaging onset time (r) itself, with 
the error being larger at early times. This is demonstrated in Figures 2D and 2E for 
the SDD and MRG models respectively, using St/t = 10%. Furthermore, we see in 
Figures 2C-E, that (jSJ) agrees very well with our numerical results. 

We can draw three conclusions from the above results. First, pre-steady-state 
morphogen interpretation inevitably leads to an increase in measurement error due to 
TAW fluctuations. Given the high precision with which gene expression boundaries are 
typically defined [26] , TAW fluctuations could be a significant impediment to pre-steady- 
state measurement. Note that this source of error is completely absent in steady-state 
where gi = for both models. Second, if a system has significant TAW fluctuations, the 
MRG model provides more precise positional information than the SDD model. Third, 
systems that interpret morphogen gradients in pre-steady-state must have additional 
regulatory mechanisms to ensure that fluctuations in the averaging onset time are well- 
controlled. 

3.3. External Fluctuations 

We now focus on the positional error due to external morphogen production fluctuations. 
The positional error due to such fluctuations, u ext , has been studied previously [23J and 
here we outline the results relevant to our conclusions. The positional error is (see 
Appendix B for details) 

ul xt « hiXi(SJ/J) (6) 

where h sdd = [1 + Xd x f(x T ,r)(l - f(x T , r))" 1 ]- 1 and h mrg = 1 (and 5 J/ J = Sj /j 
for MRG model). In steady-state h s dd = 1 and the two models have equal error. 
As with TAW fluctuations, u™ t 9 is independent of the threshold position whereas 
u^f is sensitive to the threshold position and the kinetic parameters. Since uj™t is 
independent of position, MRG-like models could be favoured in pre-steady-state when 
multiple threshold boundaries are required. 

We first investigate numerically how the positional error due to external production 
variations depends on when measurement occurs (see Appendix A). In Figure 2D we 
demonstrate that, for the SDD model with 5 J / J = 0.1, the positional error is reduced by 
pre-steady-state measurement in agreement with [23] . This increased robustness to such 
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fluctuations is a major advantage of pre-steady-state measurement in the SDD model. 
In comparison, Figure 2E shows that the MRG model (with 5 jo /jo = 0.1) has constant 
positional error regardless of when measurement occurs. Moreover, the SDD model is 
typically more robust than the MRG model to source fluctuations, regardless of the 
size of the morphogen production fluctuations, since h s dd < h mrg (with X s dd = A mr9 ), 
see Figure 2F. Importantly, in all cases, we see excellent agreement between ([6]) and 
Figures 2D-F. Clearly, embryos with very large external production fluctuations will 
favour diffusive morphogen propagation with pre-steady-state measurement. 

Morphogen gradients described by the SDD model benefit from pre-steady-state 
decoding when there are large external morphogen production fluctuations. However, 
such systems have increased sensitivity to TAW variations. Therefore, even in systems 
without significant internal fluctuations, it may still not be advantageous to interpret 
morphogens prior to steady-state. Furthermore, in systems with considerable internal 
fluctuations, the SDD model can be particularly imprecise, discouraging pre-steady-state 
morphogen interpretation regardless of the external fluctuations. The MRG model is 
not favoured by morphogens with large external fluctuations. However, the MRG model 
displays increased robustness to TAW fluctuations and it is less sensitive to internal 
fluctuations at early times. If a developmental system has to interpret morphogens 
that are direct transcription factors prior to steady-state (the embryo needs a position 
quickly but not precisely) then an MRG-like model will likely be more effective than an 
SDD model so long as external fluctuations are not the dominant source of error. 

4. Experimental Application 

Having built-up a clear understanding of the precision of pre-steady-state morphogen 
gradients in the presence of fluctuations, we can now relate our theoretical approach to 
experiment. We consider two general developmental systems: first, when the morphogen 
is a direct transcription factor, such as Bed; and second, morphogens interpreted by 
cell-surface receptors, such as Nodal. We account for TAW variations and internal and 
external fluctuations simultaneously to calculate the total positional error, e (details in 
Appendix A). The precision given by e is a lower bound on the error in determining 
xt- We have neglected other sources of error, such as noise from transcription and 
translation of target genes [42J and cell-to-cell variability [13]. Furthermore, we assume 
the measurement process is perfect. Recent experiments on the Bed-hunchback (hb) 
signaling pathway in Drosophila suggests that diffusive input noise dominates over 
transcriptional/translational noise [12], supporting our approach. 

4-1. Morphogen Transcription Factors 

Morphogens that are transcription factors are measured at nanometre-scale DNA 
binding sites (Ax ~ 3nm), so, as mentioned above, internal fluctuations can be 
significant sources of error [281 126] . We focus on the Bcd-hb signaling pathway in 
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Drosophila where hb is a target gene of Bed [231 ESI ESI- This network provides precise 
positional information, with the hb expression domain boundary defined to within 2% 
of embryo length (EL) (25j [28j [291 EE]. Using the known kinetic parameters for Bed 
[35| [26] , early measurement results in large positional error, primarily due to internal 
fluctuations, Figure 3A. This suggests that early Bed interpretation (as soon as cycle 9, 
as proposed by [231 [33] ) could well be disfavoured. The experimentally observed diffusion 
constant for Bed (D = 0.3/xm 2 s _1 [35]) and the predicted value from Bergmann et. al 
(D = l/im 2 s _1 [23J) are both insufficient to allow accurate positioning at xt ~ 240/im 
after 90 minutes, see Figure 3B (for each D, J and fi are chosen to leave the steady- 
state Bed profile unchanged). Our results include spatial averaging [26] (see [HI H5] 
for possible biochemical mechanisms for spatial averaging), without which internal 
fluctuations further penalise pre-steady-state decoding. In conclusion, Bed does not 
appear well-suited to early interpretation. 

However, Bed does not independently determine the threshold position of its target 
genes. In the early Drosophila embryo a gap-gene network regulates positioning |46j. 
We analysed a simplified model of the system (see [231 E2] and Appendix C) to check 
that our conclusions for the SDD model hold for more complex networks. In this 
model, a pre-steady-state Bed gradient defines the positions of the gap genes, and 
these spatial positions are stabilised by gap gene interactions [23]. Importantly, once a 
threshold position is defined, it does not alter significantly even though the morphogen 
concentration continues to increase. In Figure 3C we show the Hb concentration profiles 
from 20 simulations with fluctuations in J only (red lines) and all three considered types 
of fluctuations (green lines). Figure 3C clearly shows that only considering external 
morphogen production fluctuations (with corresponding positional error e s dd ~ 2%EL) 
ignores significant sources of positional error {e S( id ~ 6%EL including all fluctuations). 
Hence, even in this more complex model it seems unlikely that Bed can provide 
accurate positional information after 90 minutes (the timescale for hb expression) with 
D = lfim 2 s~ 1 . 

The above conclusions are deduced from a diffusive model of Bed morphogen 
gradient formation. Recent experiments suggest that the bed mRNA may not be 
localised in the anterior region of the embryo, but instead the bed mRNA forms its 
own gradient across the embryo [El [31]. The MRG model is a representation of such 
spatially distributed morphogen production. From Figure 3A, we find that for the MRG 
model pre-steady-state measurement does not grant improved performance over waiting 
for steady-state. In particular, the SDD and MRG models have similar positional error 
after 90 minutes. Therefore, we predict that even if Bed is driven by an underlying 
mRNA gradient, pre-steady-state interpretation is not advantageous. The MRG model 
could perform better by increasing the rate of degradation, but the parameters are 
constrained by the known timescales for gradient formation and the Bed steady-state 
profile. The MRG model can be used to describe Bed gradient formation within the 
gap-gene system described above. Again, we find that the error in positioning of the 
hb expression domain boundary is too large when interpreting the MRG model prior 
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to reaching steady-state (positional error of e mrg ~ 3%EL when only j varies but 
e mrg ~ 6%EL when we consider all three types of fluctuations). 

4-2. Morphogen Interpretation at a Cellular Level 

A second class of morphogens we consider are signaling molecules initially sampled at 
the cellular level by cell-surface receptors. In this case, there are two general stages to 
morphogen interpretation. First, the morphogen gradient must be interpreted precisely 
by cell surface receptors. Second, the receptor signal must be precisely decoded into a 
corresponding target gene response within the cell's nucleus. Typically, for morphogens 
that are members of the transforming growth factor-/? family (such as Nodal and Activin 
[171 HE]), the binding of morphogen to a cell-receptor results in the phosphorylation of 
a member of the Smad-family of transcription factors [17] 08] . Stable phosphorylation 
of the appropriate Smad protein is achieved within 15-30 minutes of exposure to the 
morphogen-bound receptor [IT]. However, the corresponding timescale for target gene 
expression (for example, Xbrachury [20]) is typically about an hour. 

Activin passively diffuses through the embryo over a scale of 300/im in a period of 5 
hours [36] (D « 2.5fj,m 2 s~ 1 ) and its concentration profile in steady-state is exponentially 
decaying [37]; the SDD model should therefore provide a good picture. From the above 
discussion, we see that interpretation of Activin by cell surface receptors occurs relatively 
quickly, therefore the condition r m <C r should approximately hold. From Figure 3D, 
we see that precise pre-steady-state decoding is possible in such a scenario due to the 
effects of internal fluctuations on Activin interpretation being small (since the effective 
cell- wide measuring volume is large). The predicted precision is similar to that found 
in steady-state. Note that the interpretation of Activin is not an on-off reading (as 
occurs in the Bcd-hb system and as assumed by our calculation). Rather one gradient is 
transformed into another. Nevertheless our qualitative conclusions, namely that internal 
fluctuations are unimportant, will be unaffected by this distinction. Recent experiments 
on the Nodal signaling pathway in zebrafish have also demonstrated clear pre-steady- 
state behaviour, with the Nodal-induced gene expression of no-tail spreading with time 
to over 200/xm in length [22]. Since Nodal is interpreted by cell-surface receptors, we 
predict that the effects of internal fluctuations are unlikely to be a significant source 
of error in no-tail boundary precision. Finally, Figure 3E demonstrates that the MRG 
model can also provide precise positional information for morphogens interpreted at a 
cellular level. 

We now consider the downstream expression of target genes. Using (0J) and typical 
parameters estimated from the Bed transcription factor, we find that an averaging time 
of about an hour is required for the downstream transcription factors to reduce the 
positional error to less than the width of the cell itself (assuming there is no spatial 
averaging). This calculation is consistent with the hour timescale for target gene 
expression. 

We conclude that pre-steady-state decoding of morphogens interpreted by cell- 
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surface receptors can occur. The difficulties posed by large internal fluctuations can 
be avoided by (i) sampling the morphogen at a cellular level, thereby ensuring a large 
measurement volume and (ii) using a long averaging time to reduce fluctuations in the 
downstream transcription factor. 

5. Discussion and Conclusions 

We have calculated the precision of positional information specified by pre-steady-state 
morphogen gradients due to internal and external fluctuations and variations in the onset 
time of measurement. Our results allow us to make the following general conclusions. 
First, morphogens that are direct transcription factors, and which are interpreted prior 
to steady-state, face considerable drawbacks due to large internal fluctuations. Second, 
even in systems without significant internal fluctuations (such as morphogens interpreted 
by cell-surface receptors), pre-steady-state interpretation is not necessarily favoured. 
The properties that confer increased robustness to external morphogen production 
fluctuations (small D, early measurement) are precisely the properties that increase 
sensitivity to TAW fluctuations (since the system is further from steady-state). So, 
when does it pay for a developmental system to 'rush' morphogen interpretation? We 
disagree with Bergmann et al. [23] that systems interpret morphogens in pre-steady- 
state to increase robustness to fluctuations. Rather, pre-steady-state precision is either 
reduced or is roughly the same as in steady-state. Therefore, alternative explanations 
for pre-steady-state decoding, such as the benefit of defining multiple genes at similar 
spatial positions [20] or the need for rapid (but potentially imprecise) gene boundaries 
are more consistent with our analysis. 

In this study, we focused on three-dimensional systems. In this case, larger detector 
regions do reduce the effects of internal noise (see Ref. [IS]). Conversely, strictly in 
two- dimensions the size of the detector region has little effect on positional precision 
[T5] . For morphogen interpretation to occur effectively in two-dimensions, the depth 
of the system needs to be comparable to the size of the morphogen protein (i.e. a 
few nanometres). Even in relatively flat systems such as eye and wing imaginal disc 
patterning in Drosophila [IHl [H3 EH [T7] , morphogen interpretation therefore effectively 
occurs in three-dimensions. Indeed, experiments on flattened Xenopus animal caps 
(only a few cells deep), find no appreciable effects on the precision of the Activin 
gradient [52]. However, it would be interesting to investigate the positional precision of 
morphogen gradients constrained to membranes as two-dimensional effects may then be 
more pronounced. 

In our analysis, we also assumed that the morphogen gradient was sampled 
continuously in space across the system. However, in real systems the morphogen 
gradient is sampled discretely at nuclei or cell positions. Moreover, in the Bcd-hb system 
in Drosophila, the number of nuclei changes during the relevant developmental period, 
doubling at the end of each cycle. For this reason, it is revealing to scale the positional 
error by the internuclear distance, to see, for example, whether neighbouring nuclei 
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can be distinguished. Since the internuclear distance during cycle 13-14 is around 
2 — 3%EL, we see from Figure [3]A. that in pre-steady-state it is never possible to 
distinguish neighbouring nuclei at mid-embryo using the Bed gradient. 

We have applied our results to several experimental systems. In particular, we 
focused on the Bcd-hb system in Drosophila which has provoked controversy regarding 
whether Bed is interpreted prior to steady-state [331 El] ■ If Bed propagation is described 
by the SDD model (as proposed by [231 [33] ) , we find that pre-steady-state interpretation 
results in increased positional error given the measured Bed diffusion constant, even 
when we include the effects of spatial averaging. We have also seen that our conclusions 
hold when we include more complex dynamics, such as the gap-gene regulatory network. 
Similar results hold for the MRG model, with early measurement of Bed disadvantaged 
by internal fluctuations. Hence, we predict that Bicoid is unlikely to be interpreted 
prior to reaching steady-state. Activin and Nodal are, however, both consistent with 
our predictions for morphogens that can be interpreted prior to reaching steady-state. 
Since the relative size of a cell is large, cells can interpret morphogen gradients quickly 
as the effects of internal fluctuations are small. Moreover, the cell then uses longer 
averaging times for the downstream transcription factors to reduce the effects of internal 
fluctuations on target gene expression. Finally, we predict that systems employing pre- 
steady-state measurement may have additional regulatory mechanisms to ensure that 
fluctuations in the onset of averaging are reduced. 

How morphogen gradients provide such precise spatial information has long 
been puzzling. Whilst pre-steady-state decoding of morphogen gradients does confer 
increased robustness to some fluctuations, it is certainly not a 'cost-free' solution for 
developing embryos, particularly in the presence of significant internal fluctuations. 
Furthermore, we have demonstrated that fluctuations in the onset time of averaging can 
have as large an effect on precision as other sources of noise. By providing a quantitative 
approach to calculating the positional error due to multiple sources of fluctuations, we 
have outlined a framework for future experimental and theoretical studies of pre-steady- 
state morphogen gradients. 
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Appendix A. Calculating Positional Error 

We find the concentration over a single realisation numerically, 



where p T , rm denotes the average morphogen concentration profile over a period r m 
starting at time r after initiation, with t n = 5t + t n -i {t\ = r), 5t = T in d and N = r m /St. 





n=l 
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For a particular gene expression boundary xt, we first find the corresponding threshold 
concentration (in the absence of fluctuations) px = Pr,T m i. x T)- For each simulation J, j 
and r are drawn from Gaussian distributions with standard deviations S J, Sjo and St 
respectively. (1A.11) is then used to calculate the concentration profile (without internal 
fluctuations) averaged over a period r m for each parameter realisation, with the error 
determined from finding where this profile goes through p?. The positional error due 
to TAW and external fluctuations is found from 50 repetitions of the simulations to 
build the distribution that describes the error in determining xt- In Figure 2, we only 
investigate TAW fluctuations or source fluctuations in each subplot (so, for example, in 
Figure 2C SJ = Sj = 0). 

To calculate the positional error due to internal fluctuations for each given 
parameter realisation we use (j4j). The constant A = k T j AxD N spat , where N spat 
denotes possible spatial averaging [26J and k T is a system-dependent constant associated 
with time averaging (see [15] for details). When calculating the total positional error, 
we simultaneously calculate the positional error from internal fluctuations u int using (Jlj) 
and the error from TAW and source fluctuations, u>TAw,ext- We find the total positional 
error due to internal noise together with TAW and morphogen production fluctuations 
using e ~ (^TAWext ^^int) 1 ^ 2 (approximation valid if T m <C r). For all results presented, 
the mean measurement error and standard deviation were found from 10 independent 
computations of the positional error. When is very small, N becomes very large 
and we approximate the sums in (jlj) and (lA.lj) by integrals. 

Appendix B. Derivation of ujtaw an d w ex t 

ujtaw'- For small fluctuations in r, we can Taylor expand the concentration at xt to 
first order with respect to fluctuations in time, Sp ~ \dtp\ST where Sp is the fluctuation 
away from the mean concentration value at position x and time r. To leading order, 
uj ~ \l/d p x\Sp pj)] and hence ujtaw ~ \dtp/d x p\ST. Substituting in (j2j) and (j3j) we 
obtain (jSj) and g s dd and Qmrg- ^ext- Similar to above, we Taylor expand the morphogen 
concentration with respect to fluctuations in J (for SDD model) or j (for MRG model) 
to find u ext w \l/d p x\Sp w \djp/d x p\SJ. 

Appendix C. Gap-gene model 

The results presented in Figure 3C are based on the simplified reaction-diffusion model 
of the Drosophila gap gene regulatory network described in [23]. The model incorporates 
the gap genes giant, knirps and Kruppel and their mutual repression of each other and 
hb [321 [46] . The Bed morphogen gradient evolves with time, allowing an investigation 
of how the system responds to a pre-steady-state morphogen gradient. The gap genes 
and hb are activated by the Bed gradient at different threshold concentrations. The 
gap genes are activated, on average, 90 minutes after Bed is initiated, consistent with 
experiment [23J. Repression, activation and autoactivation (gap genes promoting their 
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own expression) are modelled using Hill functions. The diffusion rates for the products 
of the gap genes and hb are significantly less than the Bed diffusion rate (i.e. they are 
effectively localised) but their production and degradation rates are significantly larger 
than the respective rates for Bed. Full details of equations and parameters used can be 
found in [23], though we use X s dd = 100/xm. In the simulations we alter r by fluctuating 
the gap gene activation time. We also alter r m , but since r m <C r such fluctuations 
are not a major contributor to positional error. External fluctuations are incorporated 
by fluctuating the Bed production rate, J. For each set of parameters, the threshold 
position is defined as where the normalised Hb concentration equals 0.5 after time r+r m . 

Whilst incorporating fluctuations in J, r and r m is straightforward, dealing with 
stochastic fluctuations in the concentrations is more difficult. One solution would be to 
solve the set of differential equations stochastically, a complicated procedure. Instead, 

— 1 /2 1 /2 

we fluctuate the threshold concentrations for each protein with deviation Ar m p T 
and p(xt, t) = pr- Effectively, instead of having fixed threshold concentrations and a 
fluctuating morphogen gradient, we have a smooth morphogen gradient but fluctuating 
threshold concentrations. To demonstrate that this approach is reasonable we consider 
the error due to internal noise in steady-state using the SDD model. In the former case, 
using (jlj), uj int « Ar-m 1 ^ 2 (p(xt)} 1 ^ 2 /\(d x p(xT)}\ = AX/ ^/r m p T whereas for the latter 
scenario, u int ~ \5px/pT (for 5pr -C px)- Therefore, if 5px = Ar m p T <C Pt the 
two approaches give the same positional error. This result extends straightforwardly to 
pre-steady-state measurement. 
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Figure 1. Pre-Steady-State Profiles and Fluctuations. A: Average morphogen 
concentration (scaled by steady-state value) at fixed position xt against time after 
initiation (on log scale). Black line corresponds to the MRG model whilst dashed red 
line denotes the SDD model, po = pi = 30/im~ 3 , D = lpm 2 s~ 1 , p = v = 10~ 4 s _1 , 
xt = 240/Ltm, X s dd = Xmrg = 100 pm. Inset shows ratio of the concentrations in the 
SDD and MRG models at x = xt for varying times after initiation. B: Particle number 
variance, a 2 , from simulations of SDD model (o) with line denoting a 2 — (n(x)), 
the mean particle number, the prediction from Poisson statistics. Inset is same but 
for MRG model. \ s dd = A mrg = 2.6pm, J = lpm~ 2 s~ 1 , jo = 0.67 pm~ 3 s~ 1 , 
D = Dq = 0.67/im 2 s _1 , p = v = 0.1s , r = 5s, r m = Is, Ax = \Q~ 2 pm. 
C Positional error u>i nt from stochastic simulations [15] compared with theoretical 
prediction for SDD model (HJ. Results for range of parameters: 5 < r < 60 seconds, 
1 < Tm < 20 seconds, 2 < x T < 7pm, D = D = 0.67^m 2 s _1 , 0.01 < p < 0.1s" 1 , 
1 < J < 10pm~ 2 s~ 1 . All simulations are three-dimensional, with L x — 20pm and 
L y = L z = 3pm. 
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Figure 2. Positional error due to fluctuations. (MRG Model: black lines and 
symbols. SDD Model: red lines and symbols). A: Positional error due to internal 
fluctuations, LUi nt , against averaging initiation time r. Solid, dashed and dotted lines 
correspond to D = 1, 1.5, 2pm 2 s~ 1 for the SDD model and v = 1, 1.5, 2 x 10~ 4 s _1 for 
the MRG model respectively. J, jo and p, are also altered such that each simulation 
has the same steady-state profile with po = pi = 30/ito~ 3 , A = 100/im. B: Morphogen 
concentrations for SDD and MRG models plotted at 30 minute intervals (r = 30, 60, 
90 minutes) with a fixed threshold concentration 2p,m~ 3 (dotted line). The arrows 
denote how the measured threshold position varies with time in the SDD and MRG 
models. C: Positional error due to fluctuations in r, lutaw, as a function of St/t. 
D: Comparison of positional error (°) and cu^f (filled squares) against 2Dt/x\ 

with \ s dd constant [x\/2D is the approximate time for the system to be in steady-state 
at position xt)- E: Comparison of positional errors (°) an< ^ ^""ext (filled squares) 

against vt. For vt > 2.3, the concentration is greater than 90% of the steady-state 
concentration. F: Positional error due to external morphogen production fluctuations 
uj ex t as function of the relative magnitude of the source fluctuations 6 J / J (for SDD 
model) and 5 jo/ jo (for MRG model). In Figures C-F, solid lines are theoretical 
fits from ((5]) and (J6j) . Parameters used, unless stated otherwise, xt — 240/im 
[25] . t — 85 minutes [23], r m = 5 minutes, L — 500p,m [25], J = 0.3p,m~ 2 s~ 1 , 
jo = 3 x 10 _3 /im _3 s _1 , D = lp,m 2 s~ 1 [23], p, = v = 10 _4 s _1 , X s dd = A mrg = lOO^rn 
[25], k T = 0.56 PU, Ax = 3x10" 3 ^to [26], D = 0.3pm 2 s- 1 [35], N spat = 0.06r m [26], 
SJ/J = 0.1, Sjo/jo = 0.1, St/t = 0.1. 
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Figure 3. Total positional error. (Same parameters and notation as Figure 2 
unless stated otherwise). A: Total positional error scaled by system length e/L (solid 
lines) against averaging initiation time, r. Dashed lines correspond to the positional 
error solely due to internal fluctuations (TAW and external fluctuations not shown for 
clarity). B: Total positional error in SDD model scaled by system length e s dd/L against 
rate of diffusion D with r = 90 minutes, keeping A and po = JX/ D constant so that 
we compare with the same Bed steady-state profile (line denotes the steady-state value 
of e). C: (Normalised) Hb concentration against relative embryo position in simplified 
gap-gene regulatory network model (see Appendix C and [53] for details). Red lines 
correspond to fluctuations only in J (5 J/ J = 0.1). Green lines include all three sources 
of error considered (parameters for Bed as A). Arrows denote positional error when 
relative Hb concentration is 0.5. Inset: concentration profiles of Hb (yellow) and gap 
gene products [32] (red=knirps, green=kruppel, hlue=giant). D: Same as A except 
Ax — 10/im, D = Dq = 2.5/j,m 2 s" 1 , xt — 280/xm [36j and no spatial averaging, with ■ 
and o denoting the positional error due to external and TAW fluctuations respectively 
in the SDD model. E: Same as D except for MRG model with v — 2.5 x 10~ 4 s -1 . 



